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We present calculations of the electronic and thermal transport properties of graphene antidot 
lattices with a finite length along the transport direction. The calculations are based on the 7r-tight- 
binding model and the Brenner potential. We show that both electronic and thermal transport 
properties converge fast toward the bulk limit with increasing length of the lattice: only a few 
repetitions (~ 6) of the fundamental unit cell are required to recover the electronic band gap of 
the infinite lattice as a transport gap for the finite lattice. We investigate how different antidot 
shapes and sizes affect the thermoelectric properties. The resulting thermoelectric figure of merit, 
ZT, can exceed 0.25, and it is highly sensitive to the atomic arrangement of the antidot edges. 
Specifically, hexagonal holes with pure zigzag edges lead to an order-of-magnitude smaller ZT as 
compared to pure armchair edges. We explain this behavior as a consequence of the localization of 
states, which predominantly occurs for zigzag edges, and of an increased splitting of the electronic 
minibands, which reduces the power factor S 2 G e (S is the Seebeck coefficient and G e is the electric 
conductance) . 

PACS numbers: 85.80.Fi,73.23.Ad,68.70.Lm,68.65.Cd,73.21.Cd,63.22.Rc,72.80.Vp 
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I. INTRODUCTION 

Ideal thermoelectric materials conduct electricity very 
well while the heat conduction is poor. Their appli- 
cations include power generation and refrigeration^. 
The optimization of thermoelectric properties has been 
a topic with wide interest! 3 -^, and in particular nanos- 
tructured materials for thermoelectrics is a rapidly ex- 
panding field of research. One proposal has been to 
increase the Seebeck coefficient, 5, by reducing the di- 
mensionality of the systenpHZJ Another idea is to uti- 
lize the low thermal conductance, together with a sharp 
resonance in the electronic conductance G ei of molecu- 
lar junctions^UJ Reduction of the thermal conductiv- 
ity in nanostructured materials may be achieved using 
nanomesh structures, surface-disorder and -decoration, 
passivation, or by other means. However, the electronic 
conductance should ideally not be affected. Examples 
of nanostructured thermoelectric materials include pas- 
sivated Si nanowires^, Si antidot latticeff^S, chevron- 
type graphene nanoribbons^, and connected capped car- 
bon nanotubes^H Here, we turn the attention towards 
graphene antidot lattices (GALs), a nanomesh of holes 
in graphene with promising electronic properties such as 
a tunable band gapM 

The efficiency in converting temperature gradients into 
an electric voltage, at an average temperature T, is 
quantified by the dimensionless figure of merit ZT = 
S 2 G e T ' J 'ft, where high ZT implies a good thermoelectric. 
We thus seek a high electronic power factor, S 2 G ei and 
minimal thermal conductance k = K p h + K e which in- 
cludes contributions both from phonons and electrons. 
Thermoelectric materials with ZT « 1 have an effi- 
ciency in the range of available thermoelectric compo- 
nents based on nanostructered bulk materials, whereas 
ZT > 3 is needed to compete with conventional refriger- 



ators and generator J 18 | 19 l 

Graphene can sustain current densities six orders of 
magnitude larger than copper, has a measured record 
high stiffness, and is foreseen to have numerous applica- 
tions ranging from nanoelectronics, spintronics and na- 
noelectromechanical devices 2 ^. Graphene is furthermore 
one of the best thermal conductors knowrPffiU. It has 
been predicted to posses a giant Seebeck coefficient when 
gated by a sequence of metal electrodes^. However, ways 
to reduce the superior thermal conductivity of graphene 
are needed if one looks for thermoelectric applications. 
Several ways to reduce the thermal conductivity have 
already been examined, such as interface mismatching 
betwee n grap hene and nanorib bon d 2 4 * 25 !, the presence of 
isotopes^Mm anc [ point defect <PQ|31J Edge disorder has 
been predicted theoretically to suppress heat conduc- 
tance of graphene nanoribbon o 32 ^ 33 ^ and ZT exceeding 3 
has been theoretically predicted for such systems in the 
diffusive limit 34 . 

Graphene antidot lattices (GALs) have been proposed 
as a flexible platform for creating a semiconducting mate- 
rial with a band gap, which can be tune d by v arying the 
antidot size, shape, or lattice symmetr ^ 17 * 35 ^ 3 -^. GALs 
can be fabricated by electro n bea m lithograph ^ 38 * 39 *, 
by block copolymer lithograph ^ 40 * 41 I with hole distances 
down to 5nm, and at a larger scale through nanorod 
photocatalysis 42 and anisotropic etching 43 . To the best 
of our knowledge, no studies have been reported on the 
thermal properties of finite GALs. Apart from their in- 
trinsic scientific interest, these studies are necessary to 
assess whether the thermal properties can be engineered 
in a manner similar to the electronic case. Of course, 
all realistic devices are of finite length, and the study 
of size-effects is important for practical purposes. For 
completeness, we mention here other related studies that 
have recently been reported. A number of studies of 



electron and/or phonon transport properties of regular 
defects in ribbons are available, see, e.g., Refs. [44H46] . 
Recently, Lopata et al. studied electron transport of in- 
finite GALs 47 . Finally, during the preparation of this 
manuscript Karamitaheri et al. reported a combined 
study of electron and phonon transport properties based 
on the band structures of infinite G ALs 4 ^ and Tretiakov 
et al. reported results for topological insulators 49 which 
share certain key properties (e.g., flat bands) with GALs. 

The topic of this paper is thus the electronic and ther- 
mal transport properties of finite graphene antidot lat- 
tices. The finite GALs are viewed as a part of an inte- 
grated graphene-based system, e.g. used as an electrode 
for molecular conductors'^, see Fig. [I] In order to shed 
light on the question to what extent it may be possible to 
engineer the thermoelectric properties we investigate how 
different antidot shapes and sizes affect the thermoelec- 
tric properties. Interestingly, even though the base ma- 
terial - graphene - is an outstanding thermal conductor, 
we find that the resulting thermoelectric figure of merit, 
ZT, can exceed 0.25. However, ZT is highly sensitive to 
the atomic arrangement at the edge of the etched holes, 
partly due to electronic quasi-localized edge-states. As 
we shall show below, this favors antidots with armchair- 
type edges for thermoelectrics. 

The paper is organized as follows. In section [IT] we 
introduce the systems and outline the theoretical and 



numerical methods used. In section III we present our 
results for the electronic transport properties of GALs. 
Especially, we discuss the convergence with number of 
repetitions of the basic unit cell, and also consider the 
localization of electronic states at zigzag edges. In Sec- 



tion IV we examine the influence of the perforation and 
geometrical effects on the thermal transport properties. 
This leads to the analysis of thermoelectric properties in 
section [Vl The results are summarized and discussed in 
section [VTl 



II. SYSTEMS AND METHODS 

Throughout in this paper we focus on triangular 
graphene antidot lattices: these syste ms ar e known to 
lead to a gap in the electronic spectru m 136 1 37 1 which is es- 
sential for the present purposes. Due to the high lattice 
symmetry the number of independent lattice parameters 
is small, and furthermore, these systems are the most 
thoroughly studied, both theoretically and experimen- 
tally. Recent experiments have illustrated that hexago- 
nal antidots may stabilize with pure zigzag and armchair 
edge chirality^. The antidot diameter, shape, position 
and the ratio of removed atoms to unit cell size are all 
important parameters which we examine to get a full pic- 
ture of the electronic, thermal and thermoelectric prop- 
erties of GALs. Another important variable is the length 
of the region exposed to the nanoperforation. The sys- 
tems studied here consist of an antidot lattice of finite 
length connected to two pristine graphene leads (Fig. [T|), 



and the infinite direction perpendicular to the transport 
direction is treated using periodic boundary conditions 
and corresponding /c-point sampling. 

M=2 
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FIG. 1. (Color online) System setup and the computational 
rectangular unit cell (green rectangle). Two graphene leads 
are connected by the finite GAL. The depicted system is a 
{10,5zz} GAL with a length of 2 (M = 2) corresponding to 4 
holes in the direction of transport. 

We use the nomenclature introduced in Ref. [17] and 
specify the GAL by {Lsi,Ss2}, where the L is the length 
of the side of a hexagonal Wigner-Seitz cell, and S is the 
side length of the antidot (see Fig. [I]), both in units of 
the lattice constant a/3&o, with ao = 1.42 A being the 
carbon-carbon distance. The label S2 = {zz, arm, cir} in- 
dicates whether the hole has zigzag or armchair edges, 
or if the hole is circular resulting in mixed armchair and 
zigzag edges. In principle, one could also consider differ- 
ent sheet orientations compared to the transport direc- 
tion; we do not present a systematic study here, and fix 
the underlying graphene sheet as armchair (s\ = arm, as 
in Fig. TJ. We have tested a selection of 'zigzag sheets' 
and did not find any qualitative difference with respect 
to thermoelectric properties. As an example {10,5zz} is 
a L = 10 antidot lattice with transport direction perpen- 
dicular to the armchair direction^ and with a hexagonal 
hole with same orientation as the lattice hexagons result- 
ing in zigzag edges and a side length of S = 5 (see Fig. pi). 
Armchair edges are obtained if the hexagonal holes are 
rotated opposite to the lattice hexagons (see also Fig. [2] 
below). 



A. Method 

Both electronic and phonon transport properties 
are calculated from atomistically determined energy- 
dependent transmission functions, T e and T v hi as de- 
scribed below, and using these in a Landauer-type for- 
mula. For a spin degenerate electronic system the Lan- 
dauer formula reads, 
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where uf(E, jjLl/r) is the Fermi-Dirac distribution at the 
chemical potential of the left /right lead. We will employ 
this in the linear-response limit, and consider variations 
with changes in the chemical potential, e.g. by doping or 
gating of the graphene system. The following integrals 
can be evaluated from the electronic transmission, 
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They relate the electronic current and the electron heat 
current, Iq, in the linear response regime: 
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where A/i = [ii,— jir and AT — T^ — Tr. From these in- 
tegrals several physical properties follow 52 ; the electrical 
conductance G e {ji) = Jy = e 2 L , the electron thermal 

conductance K e {ii) = L2 — j^ /T, and the Seebeck co- 
efficient S(fi) = ^|/ e =o = ^r- 

For phonons the Landauer formula takes an analogous 
form, 
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where ub(oj) is the Bose distribution function. Again 
we use it in linear response and consider the thermal 
conductance from phonons given by, 
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Both transmission functions, T e and T p h are obtained 
using a recursive Green's function method (see Ref. [53] 
and the references cited therein) with self-energies rep- 
resenting the semi-infinite perfect graphene electrodes. 
The self-energies, ^l,r, are iteratively constructed from 
the semi-infinite graphene left (L) and right (R) leads. 
The calculation of both electron and phonon /c-averaged 
Landauer transmissions together with the thermoelectric 
properties are pe rformed by an atomistic Green's func- 
tion method 16 * 54 *, 
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Here the retarded Green's function G r D (E, k), is obtained 
from the Hamiltonian, H, G r D = [EI - H - H r L - ££>] _1 , 

and the broadening matrices due to the electrode cou- 
pling are defined as Tl^r = i[S£ R — J2 L R ]. The param- 
eter Nk gives the number of sampled /c-points 55 . Similar 
equations hold for the phonon transmission: the Hamil- 
tonian is replaced with the dynamical matrix H — >• K, 
and the energy is replaced with EI — >> cj 2 M, uj being the 
frequency and M is the diagonal mass matrix. We first 



perform a structural relaxation and then calculate the 
Hamiltonian/dynamical matrix for three unit cells (M=3, 
6 holes) between the pristine graphene leads. The ele- 
ments corresponding to the center cell (2 holes) are then 
subsequently repeated to increase the length of the GAL. 

The electronic system is modeled by a nearest- neighbor 
7r-model (V pp7r = 2.7 eV) together with the Harrison scal- 
ing law to take into account the changes in the hopping 
matrix element due to the edge relaxation 56 . Based on 
the same method Guinea and co-workers^ have shown 
how strain in graphene can lead to a pseudo-magnetic 
field affecting the electronic properties. We find that the 
modulation of the hopping elements is of minor impor- 
tance for the present applications. In order to examine 
the effect of passivation we have performed calculations 
of the band structures with a model including two d- 
orbitals for each C-atom and an explicit model for the 
carbon- hydrogen interaction 5 -^. The qualitative features 
of the band diagram, and the edge states discussed be- 
low, depend surprisingly little on the presence of hydro- 
gen passivation. However, the effect of passivation might 
be more important in antidot lattices with localized zero 
energy states such as in the triangular antidots consid- 
ered in Ref. HU 

The dynamical matrix is computed using the empiri- 
cal Brenner interatomic potential 60 . This is done for the 
system cell by the finite difference approach after a struc- 
ture relaxation performed by the 'General Utility Lattice 
Program' (GULP/ 6 -^. Momentum conservation is impor- 
tant for low frequency modes and it is imposed after the 
finite difference calculation by adjusting the diagonal el- 
ements of the dynamical matrix. A few representative 
results of the relaxation are shown in Fig. [2j We see that 
the changes in the bond lengths compared to a$ are be- 
low 4%. Compression of bonds at the edge is followed 
by a region with elongation of bonds and the relaxation 
is confined in small regions in space. Matching of edge 
relaxations can result in longer-ranged relaxations (small 
compressions) emanating from the corners. This occurs 
mainly for zigzag edges (Fig. [2]). 

The results presented below are based on a set of elec- 
tron and phonon transport simulations of 20 configura- 
tions with varying hole size and lattice parameters. In 
addition to this set the electronic structure of another 27 
systems was studied to examine on the formation of edge 
states. Finally, a number of systems has been studied 
with either the zigzag transport direction or a rectangu- 
lar lattice of holes. No qualitative differences in thermo- 
electric properties were found for these systems and we 
shall not present these results here. 



III. ELECTRON TRANSPORT 

In this section we consider the electronic properties of 
finite G ALs. The result for a series of unit cell repetitions 
(M) is shown in Fig. [3] As can be seen the electronic 
transmission T e converges fast toward a length indepen- 
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FIG. 2. (Color online) Change of bond lengths due to the 
relaxation of the graphene antidots. Top: Coloring scheme 
for the bond lengths. The distribution of bond lengths af- 
ter relaxation is given for five different lattices. From top to 
bottom the considered lattices have a) small armchair holes 
{10,3arm}, b) large armchair holes {10,6arm}, c) small zigzag 
holes {10,3zz}, d) large zigzag holes {10,5zz} and e) mixed 
edges {10,4.7cir}. When the hole size is increased the spec- 
trum broaden and peaks occur at different bond lengths. The 
peaks occur at different positions characteristic of different 
hole types. The three bonds in an armchair are characterized 
by a red-pink-red color sequence. 



dent result. The behavior of the transmission function 
can directly be traced back to the band structure of the 
infinite GAL if one defines a transport band gap as the 
energy range where the transmission is below a certain 
small value. The band gap is in general found to converge 
to that of the infinite antidot lattice found from the band 
structure and the system behaves 'bulk-like' after only six 




FIG. 3. (Color online) Top: Convergence with length (M, 
the number of unit cells with two holes along the device) of 
transmission for a {10,3arm} antidot lattice. Bottom left: 
Zoom at the band gap for the {10,3arm} GAL. The leftmost 
vertical dashed line marks the value of the band gap obtained 
from the band structure of an infinite GAL. Bottom right: 
Transmissions close to the Fermi level for the selected GALs 
shown in Fig. [2] The transmission curves have been shifted 
by an integer to ease the comparison. 



to seven unit cell repetitions. Thus the transport band 
gap can be determined from a calculation of the disper- 
sion on an infinite GAL using a primitive unit cell due to 
the fast convergence property illustrated in this section. 
The converged values of the band gap are given in Fig. [4] 

To access the effect of relaxation on the electronic 
structure we have plotted the obtained band gaps ne- 
glecting the modulation of the hopping elements for the 
results with varying hole size for the armchair holes 
(black squares) and zigzag holes (black diamonds). The 
relaxation is found not to play a qualitative role in the 
equilibrium electronic properties of GALs within this 
model. The convergence is independent of the lattice 
parameters, and in all simulations presented hereafter 
between 8 and 10 unit cells are used. 

Besides the band gap we observe that it is possible 
to approximate the transmission versus energy as linear 
curves corresponding to a simple reduction of the pristine 
transmission, % oc \E\ (see example in the transmission 
plots in Fig. IsT) . We have calculated envelope lines ob- 
tained from a scaling of the pristine transmission with 
the width of the constriction T e fi = R e fi 7o, where 7o is 
the transmission of pristine graphene. The reduction fac- 
tor, i? e ff 5 describes the amount of pristine transmission 
that survives the lattice perturbation in terms of a regu- 
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FIG. 4. (Color online) Scaling of the electronic band gap E g 
with increasing ratio of removed atoms compared to the sim- 
ple scaling law estimated for circular holes 17 . The systems 
considered are {10,Szz} with S = 3,4,4.5,5,5.5, {10,*Sarm} 
with S = 3,4.5,5,6, and {10,Scir} with S = 3.5,4.7,5. 
Furthermore we include two sets, {L,3arm} and {L,3zz} for 
L — 6, 8, 10, 12, with fixed hole geometry. 



lar perforation. The actual reduction factor is estimated 
as the average reduction found at each energy point. We 
find that the electronic transmission is reduced more than 
what would be expected from the effective width reduc- 
tion R e ff = W/Wo- Here W is the minimal width along 
the device and Wo is the width of the pristine graphene 
sheet. For the systems considered in Fig. [2] the hole di- 
mension is varying between 1.2 nm and 2.6 nm giving an 
effective width reduction between 71% and 26%. The ac- 
tual reduction factor is descreasing linearly with hole di- 
mension from 24% to 5%. Therefore, only a minor part of 
the average transmission reduction can be ascribed to the 
narrowing of the conducting plane. The present model 
does not take special account of the band gap opening. 
One could instead ask if the peak transmission is limited 
by the effective width. The peak transmission reduction 
factor is found to be decreasing from 65% to 21%, and fits 
the effective width reduction very well for small holes. As 
the hole size increase the effective width is overestimated 
due to the triangular lattice structure of the perforation, 
and the reduction factor approaches the averaged value. 



A. Localization at zigzag edges 

There is an important difference between holes of dif- 
ferent atomic arrangements at the edges. In Fig. [5] we 
compare the electronic transmission and the band struc- 
tures for two large holes with zigzag and armchair edges, 
respectively. The figure illustrates how the transmis- 
sion can be directly traced back to the band structure 
of the GAL. It is furthermore seen how the structure 
with zigzag edges leads to an additional splitting into 
fiat minibands around the Fermi level. This feature can 
be understood in terms of localized states due to a lo- 
cal excess of atoms of one of the two sublattices in the 



graphene bipartite lattice^. The local imbalance of A 
and B type atoms at the edges leads to the corresponding 
number of defect states. In hexagonal holes with zigzag 
edges each side consists of a segment of either type A or 
B atoms. The hybridization will be small between these 
defect states which are partially separated in space. As 
illustrated at the rightmost of Fig. [5j the fiat minibands 
are highly localized at the edges. In the case of hexag- 
onal holes with armchair edges each side consists of an 
alternating sequence of A and B atoms. Therefore, these 
defect states hybridize more resulting in a larger shift 
from the Fermi level and a reduced flatness of the bands. 
As can be seen from Fig. [5] the first bands with minimal 
dispersion are mainly localized in the small zigzag corner 
region between two AB sequences. 

It is possible to quantify the degree of localization from 
the weight of the eigenstate at each atom. The local- 
ization factor, for a given eigenstate in th e site basis 
^n = [^i?^2? ■■■j'Mjv]? is nere defined asj 62 ! 63 * 
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This factor equals 1/N when the state is fully delocalized 
and all weights Ui have the same value. In the case of 
a state localized at a single site it gives 1. The inverse 
localization factor gives a measure of the number of sites 
that contribute to a given state. 

A numerical example is given in Fig. [6] For the sys- 
tem with L = 10 and zigzag edges we see that the first 
conduction states corresponding to the flat minibands 
are more localized than the following bands. For the 
armchair edge the localization parameter shows a weaker 
dependence on the state index. We conclude that the 
flat minibands of the zigzag edge are more localized than 
the corresponding states of an antidot lattice with arm- 
chair edges. Figure [6] furthermore illustrates that this 
conclusion is independent of the values considered here 
for the unit cell dimension L. The localization factor 
for the armchair edge depends more strongly on the di- 
mension of the unit cell, that is, the hole-hole distance, 
and is in general an order of magnitude lower compared 
to the zigzag edge. Therefore, the band gaps of GALs 
with armchair edges are determined by the confinement, 
as opposed to the case of zigzag edges where it is gov- 
erned by edge state formation. This is the reason why the 
band gap scales differently depending on the edge type 
of the hole (Fig. E|. For very small holes with zigzag 
edges we find that the almost dispersionless GAL mini- 
bands are positioned further into the band structure, see 
the {L,3zz} curve in Fig. [4] However, as the length of 
the edge is increased the zigzag edge-state energies are 
located directly at the band gap. It is therefore not to 
be expected that a larger hole with certainty results in 
a larger band gap. Even though this happens for holes 
with armchair edges, introduction of zigzag regions may 
suppress the band gap, which is important for electronic 
and optical applications of antidot lattices. 



Electronic band structure 




FIG. 5. (Color online) Left: Band structure for {10arm,5zz} (top) and {10arm,6arm} (bottom) antidot lattices respectively. 
Middle: Corresponding electronic transmission around the Fermi level. The energies of the states illustrated and compared 
further have been marked with circles. Right: Eigenstates 1, 3, 4 and 6 at the T-point with energies as marked in the 
corresponding band structures. The eigenstates of the {10,5zz} antidot lattice (top) are very localized at the edges. The 
eigenstates of the {10,6arm} antidot lattice (bottom) are less localized at the edges, but resemble corner states. A phase of 
zero and tv is colored blue and yellow respectively. 



We conclude that GALs with armchair edge geometry 
have a larger band gap as compared to both zigzag edge 
geometries and the predicted scaling^. Furthermore, the 
hexagonal antidots with armchair edges show a system- 
atic scaling of the band gap with hole size making this 
system preferable for electronic applications. 



IV. HEAT TRANSPORT 

We now turn to the thermal transport properties of 
finite GALs. In Fig. [7] we show the phonon transmission 
as a function of M, the number of repeated unit cells. To 
quantify the convergence with length we have calculated 
the thermal conductance at 300 K for the antidot lattices 
of Fig.[2]at different lengths. This shows that the thermal 



properties converge at a length scale similar to that of the 
electrons so also the phonons behaves 'bulk-like' after six 
to seven unit cell repetitions. In all simulations presented 
hereafter we use 8-10 unit cells. 

The thermal conductance due to phonons in pristine 
graphene at 300 K should be compared to a measured 
thermal conductivity of cr^ p « 4.5-5.5 xlO 3 W/(m K). 
Our result compares well to other theoretical calcula- 
tions, where it similarly was reported^ that the re- 
duced ballistic thermal conductance, which we find to 
be nlt/iW Q h) « 4.27 x 10 9 W/(m 2 K), is much larger 
than the experimentally extracted partially diffusive re- 
sult cr^ p /(L exp ) « 0.39 - 0.48 x 10 9 W/(m 2 K). Here 

h = 3.35 A and L exp « 11.5 /im are the graphite inter- 
layer distance and traveled distance by the phonons in the 
experiment by Balandin et a/.™ respectively, and Wo is 
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FIG. 6. (Color online) Localization factor as a function of con- 
duction state index for a selection of distances between neigh- 
boring holes. We compare two fixed hole geometries, namely 
the {L,5zz} (top) and {L,6arm} (bottom) GALs. These fig- 
ures illustrate the physical difference between localization of 
states in GALs with pure zigzag and armchair edges through 
the variation of L and thereby the distance between the holes. 
Localization of states in GALs with zigzag edges is due to edge 
state localization as opposed to GALs with armchair edges, 
where localization is a result of the confinement of the elec- 
trons. 
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FIG. 7. (Color online) Top: Convergence with length of 
phonon transmission at 300 K for a {10,3arm} antidot lat- 
tice. Bottom left: Convergence of the thermal conductance 
from phono ns with length, normalized by its value at M=8. 
Bottom right: Comparison of the thermal conductance as a 
function of temperature for the selected GALs shown in Fig. [2] 



the computational unit cell width. The main difference 
here can probably be attributed to isotopes, electron- 
phonon scattering and especially anharmonicity being 
important for long devices. 

In analogy with the electronic transmission we have 
calculated an average reduction factor for the phonon 
transmission. The transmission of the lowest acoustic 
and especially the highest optical modes is in general re- 
duced more than the remaining of the phonon spectrum. 
Similarly to the electron case, the average reduction fac- 
tor decreases linearly with the hole width for the consid- 
ered systems. The average phonon transmission reduc- 
tion factor is found to be of the same order of magnitude 
as compared to the electron transmission. Once again 
only a minor part of the transmission reduction can be 
ascribed to the reduction in effective width of the con- 
ducting plane due to the perforation. There is a tendency 
that for small hole sizes the phonons are scattered more 
than the electrons by the nanoperforation. Furthermore, 
the electronic reduction factor can be much larger at a 
specific chemical potential for small holes. For large hole 
dimensions both the electrons and phonons are scattered 
to an extent where the transmission is reduced by more 
than 80% on average for the systems considered. For the 
largest holes up to 36% of the atoms have been removed 
from the pristine graphene plane. 

In Fig. [8] the temperature and hole size dependence of 
the phonon thermal conductance is given for our selec- 
tion of systems with varying hole size and shape. Fig. [8] 
illustrates how the thermal conductance decreases al- 
most linearly with the hole size for typical perforation 
removal ratios (larger than 5% perforation). Further- 
more, the graph shows that the thermal conductance 
has a tendency to be slightly larger for holes with zigzag 
edges (shown as diamonds in Fig. M) . A similar behav- 
ior has been found for graphene nanoribbons with zigzag 
edges 45 . However, compared to the electronic case the 
thermal transport features are less sensitive to the exact 
shape and edge of the holes. 

For the purpose of making it easy to compare our re- 
sult with other calculations and experiments we give an 
empirical expression for the thermal conductance. In the 
regime where the thermal conductance is linear in the 
hole dimension one can parameterize the thermal con- 
ductance as 



tlph 



-a(T) 



N r , 



N t , 



+ P)<h(T). 



(8) 



From this approximation we fit the 'linear regime offset' 
f3 « 0.25 and the dimensionless parameter a(T), given 
in the inset of Fig. [8] The lines in Fig. [8] are illustrating 
this parametrization. 

In Fig. [9] we have illustrated the electronic contribution 
to the thermal conductance at four different tempera- 
tures. Due to the vanishing electronic density of states 
around zero chemical potential the thermal conductance 
of GALs is dominated by phonons as is also the case of 
pristine graphene. However, the electronic contribution 
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FIG. 8. (Color online) Thermal conductance from phonons as 
a function of hole dimension. The red squares, blue diamonds 
and green circles label holes with armchair, zigzag and circu- 
lar/mixed edges respectively. Four different temperatures are 
plotted for each system. From top to bottom the thermal con- 
ductance is found at a temperature [450,300,150,90] K. The 
thermal conductance at these four temperatures is for pris- 
tine graphene found to be [8.5,6.1,2.6,1.2]^. Inset: fitted 
dimensionless parameter a describing the scaling with hole 
size of the thermal conductance. 
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FIG. 9. (Color online) Electronic contribution to the ther- 
mal conduction. Top: {10,6arm} antidot lattice. Bottom: 
{10,5zz} antidot lattice. From top to bottom the curves are 
found at a temperature of [450,300,150,90] K. The correspond- 
ing phonon thermal conductance has been marked to the right 
of the plot for comparison. 



can dominate - even at room temperature - when a large 
gate bias is applied. 



THERMOELECTRIC FIGURE OF MERIT 



Next we report the thermoelectric properties of the 



considered GALs. In Fig. 10 we compare the Seebeck co- 
efficient for a GAL with armchair edges (top) and zigzag 
edges (bottom) . The Seebeck coefficient displays peak 
values 65 of the order of 0.1-1.5 mV/K, which is similar to 
what has been obtained for other carbon based nanosys- 
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FIG. 10. (Color online) Seebeck coefficient at the four differ- 
ent temperatures [450,300,150,90] K for the {10,6arm} (top) 
and {10,5zz} (bottom) GAL. Notice the different scale on the 
Seebeck coefficient. The main Seebeck peaks of the GAL with 
zigzag edges occur at the energies of the low transmitting lo- 
calized states, whereas the main contribution to the thermo- 
electric power of the GAL with armchair edges is positioned 
at energies further into the band structure. 



terns and molecular contacts^ 8 * 30 * 66 *. 

For bulk material the thermoelectric figure of merit 
is defined in terms of the electrical and thermal con- 
ductivities, cr e ,<T£, as ZT = TcrS 2 /<jf For the ballistic 
graphene systems we can write it in terms of their re- 
spective conductances by introducing a width, effective 
length and thickness, ZT = TG e S 2 /(n e + K p h)- Maxi- 
mal thermoelectric figure of merit, ZT, is obtained after 
length convergence due to the increased band gap and de- 
creased thermal conductance. The obtained ZT shown in 
Fig. [TT] as a function of chemical potential has a number 
of peaks corresponding to a large variation of the trans- 
mission with energy. The Seebeck coefficient is a mea- 
sure of these changes and their robustness to temperature 
smoothening. One important feature is that the high 
peaks in the Seebeck coefficient for the {10,5zz} lattice 
mainly occur at very low energy, where the transmission 
is low, whereas for the {10,6arm} lattice the dominat- 
ing peaks occur at higher chemical potential. Therefore, 
the peak ZT is higher for the {10,6arm} as a result of 
the higher electronic conductance at peak position as il- 
The Seebeck coefficient is highly 



lustrated in Fig. 11 



sensitive to the variations in the electronic transmission 
resulting from different hole edges, sizes and so forth. 



In Fig. 12 we collect the maximum ZT values we have 
found for a selection of GALs. It seems possible to ob- 
tain larger ZT from GALs based on hexagonal holes with 
armchair edges. This is a result of the additional split- 
ting into minibands for zigzag edges. The reason for this 



additional splitting is, as mentioned in III A the forma- 
tion of edge states at zigzag edges. As a consequence the 
Seebeck coefficient can be larger for zigzag edges. How- 
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FIG. 11. (Color online) ZT for {10,6arm} and {10,5zz} lat- 
tices at the four temperatures [450,300,150,90] K. At low tem- 
perature the pure electronic figure of merit, ZT e i, can be very 
large due to a vanishing thermal conductance from electrons 
and sharp features in the transmission spectrum. At low tem- 
perature many sharp transmission features also becomes vis- 
ible in the actual ZT. 
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FIG. 12. (Color online) ZT dependence on the ratio of re- 
moved atoms in the nanoperforation at T = 300 K. Sys- 
tems included in the figure are {10,Szz} with £=3, 4, 4. 5, 5, 5. 5, 
{10,Sarm} with 5=3,4.5,5,6 and {10,Scir} with 5=3.5,4.7,5 
(full lines) and furthermore two set of systems with a fixed 
hole {L,3arm} and {L,3zz} with L=6,8,10,12 (dotted lines). 



ever, the power factor is significantly lower due to the 
lower transmission from the isolated energy levels with 
low dispersion. There is also a weak trend that the hole 
dimension compared to the system size should be max- 
imized. By increasing the hole dimension we actually 
reduce the electronic figure of merit, denned as 



ZT el = 



S 2 G e T 



Kph 



-ZT. 



(9) 



but obtain a larger fraction of it due to a reduced phonon 
conductance. Higher ZT could possibly be obtained by 
increasing the hole dimension even further, but these sys- 
tems will be very challenging to fabricate. 



The electrons-only result, ZT e \ (k^ = 0) describes an 
upper bound of the figure of merit. However, we find 
it to be somewhat artificial, due to the fact that the 
phonon contribution to the thermal conductance shifts 
the position of the peaks and ZT e \ posses a (in principle 
unbound) peak every time the electronic thermal conduc- 
tance is zero (Eq. (|9|). Especially in the presence of gaps 
in the electronic band structure, the computation of ZT e \ 
can be numerically challenging. However, evaluating the 
ZT e i expression at the true peak position can give an es- 
timate of the gain by a further reduction of the phonon 
conductance. For the {10,6arm} GAL the first peak 
(/i = 0.37eV) and the highest peak (fi = 1.23 eV) have 
a ZT = 0.17 and 0.26 with corresponding ZT ei = 4.78 
and 0.77, a factor of 28 and 3 larger than the true ZT, 
respectively. For the {10,5zz} lattice we have the first 
peak value ZT = 0.13 with corresponding ZT e \ — 5.77 a 
factor of 44 larger at the same energy. At high chemical 
potential the main limitation is the electronic structure 
and not a further reduction of the phonon heat conduc- 
tance. On the contrary one could obtain a significant 
ZT enhancement at low chemical potential by further re- 
ducing the thermal conductance. Isotope scattering, an- 
harmonic interactions, electron-phonon interactions and 
graphene-substrate interactions could all contribute to a 
reduction of the phonon thermal conductance. In Fig . 
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FIG. 13. (Color online) ZT variation with a decrease of the 
phonon thermal conductance at T — 300 K. A describes the 
fraction of the original thermal conductance used in the cal- 
culation. Top: ZT as a function of the chemical potential and 
the phonon thermal conductance for the {10,6arm} (left) and 
{10,5zz} (right) GAL. Bottom: Peak ZT as a function of the 
variation of the phonon thermal conductance. 
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13 we illustrate the effect of a reduction of the phonon 
thermal conductance. The parameter A gives the frac- 
tion of the original phonon thermal conductance kept in 
the calculation. For the {10,6arm} GAL the first peak 
(jjl ~ 0.37 eV) increases more rapidly than the high en- 
ergy peaks. The peak position crossover happens around 
A ^0.35. When the phonon thermal conductance domi- 
nates, the figure of merit variation goes as ZT/A, clearly 
present in the low energy ZT variation (see Fig. [l3j bot- 
tom) even when the phonon thermal conductance is re- 
duced below 5% of its original value. 

We propose one of two routes to obtaining higher ZT. 
Either one could find a way to reduce the thermal conduc- 
tance without affecting the electrons. Surface decoration 
might be a promising way to obtain this. Another route 
could be to improve the electronic properties of GAL's, 
e.g. by combining this system with other nanostructured 
devices. This could increase the peak ZT obtained at 
high chemical potential. We note that one-dimensional 
nanosystems may display a large Seebeck coefficient lead- 
ing to larger ZT. However, on the other hand the two- 
dimensional GAL structure may have an important ad- 
vantage. In a 2D GAL the system geometry does not pose 
inherent limits for the converted power, as it does for a 
system of parallel quasi-one-dimensional systems, such as 
quantum wires or graphene nanoribbons: the power can 
be increased by simply making the GAL-system broader. 
This provides a strong motivation for further investiga- 
tions of GALs as thermoelectric devices, perhaps as an in- 
tegrated element in future graphene nanoelectronics. De- 
spite the high intrinsic thermal conductance of graphene 
it is noticeable that one can utilize the nanoperforation 
to obtain ZT exceeding 0.25, a factor 35 enhancement 
compared to what is found for bulk graphene. 



VI. DISCUSSION AND CONCLUSION 

We have theoretically shown that GALs allow the si- 
multaneous manipulation of both electronic and thermal 
transport properties of graphene sheets. Our calculations 
have been carried out in the ballistic limit, which gives a 
reasonable first estimate for short devices whose dimen- 
sions are smaller than the various scattering lengths (im- 
portant scattering mechanisms include the anharmonic 
phonon-phonon interactions, electron-phonon scattering, 
and electron-electron scattering). Also spin-polarization 
may turn important: recent studies have shown that 
one can have spin-splitting and a magnetic moment in 



triangular^HUl antidots with pure zigzag edges. Above 
all, the most important future task is a systematic study 
of disorder effects. Our preliminary results suggest that a 
low degree of disorder can increase ZT due to a decrease 
of the thermal conductance, whereas a high degree of 
disorder affects both electrons and phonons so that the 
decrease in power factor outweights the decrease in ther- 
mal conductance. 

A key result of our analysis is the convergence of trans- 
port properties with length for GALs. The ballistic trans- 
port properties converge fast towards that of the infinite 
antidot lattice. We have also found that the quantiza- 
tion is an important feature of both electron and phonon 
transport properties of GALs. This is seen from the fact 
that the transmissions are reduced far more than what 
would be expected from an effective width estimation and 
therefore the exact scattering rate for the different edge 
types is important. The average transmission reduction 
factor is found to be on the same order of magnitude for 
electrons and phonons. In general, the formation of edge 
states determine the band gap of GALs with pure zigzag 
edges as opposed to pure armchair edges, where the band 
gap is determined by the confinement of electrons. Fur- 
thermore, the different edge characteristics play an im- 
portant role in the observed difference in thermoelectric 
properties. ZT is found to be lower for GALs with zigzag 
edges due to the additional splitting into minibands for 
large structures and a corresponding lower power factor. 
The maximal thermoelectric efficiency ZT ^0.3 has been 
obtained for GALs with pure armchair edges. Therefore, 
it is possible to obtain fair thermoelectric properties of 
graphene-based nanosystems, even despite of lattice dis- 
tortions which highly affect both the 7r-electron deter- 
mined electronic properties and the sp2-bonding deter- 
mined thermal conductance, such as the nanoperfora- 
tions. The main limitation in thermoelectric applications 
of GALs at high chemical potential is set by the electronic 
structure because the electronic heat conductance is large 
at the high energy peak position of S and ZT. At low 
chemical potential we expect that one could benefit from 
a further reduction of the phononic thermal conductance 
due to isotope scattering and anharmonic interactions. 
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